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Abstract 

We consider a general d-dimensional quantum system of non-inter- 
acting particles, with suitable statistics, in a very large (formally infi- 
nite) container. We prove that, in equilibrium, the fluctuations in the 
density of particles in a subdomain A of the container are described by 
a large deviation function related to the pressure of the system. That 
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is, untypical densities occur with a probability exponentially small in 
the volume of A, with the coefficient in the exponent given by the ap- 
propriate thermodynamic potential. Furthermore, small fluctuations 
satisfy the central limit theorem. 



1 Introduction 

Statistical mechanics is the bridge between the microscopic world of atoms 
and the macroscopic world of bulk matter. In particular it provides a pre- 
scription for obtaining macroscopic properties of systems in thermal equilib- 
rium from a knowledge of the microscopic Hamiltonian. This prescription 
becomes mathematically precise and elegant in the limit in which the size 
of the system becomes very large on the microscopic scale (but not large 
enough for gravitational interactions between the particles to be relevant). 
Formally this corresponds to considering neutral or charged particles with 
effective translation invariant interactions inside a container and taking the 
infinite- volume or thermodynamic limit (TL). This is the limit in which the 
volume |V^| of the container V grows to infinity along some specified regular 
sequences of domains, say cubes or balls, while the particle and energy den- 
sity approach some finite limiting value |^ 0, |G|, ^]. This limit provides a 
precise way for eliminating "finite size" effects. 

It is then an important result (a theorem, under suitable assumptions) of 
rigorous statistical mechanics that the bulk properties of a physical system, 
computed from the thermodynamic potentials via any of the commonly used 
Gibbs ensembles (microcanonical, canonical, grand canonical, etc.), have well 
defined "equivalent" TL's [0, 0, |3- These free energy densities are further- 
more proven to be the same for a suitable class of "boundary conditions" 
(b.c), describing the interaction of the system with the walls and the "out- 
side" of its container. When this independence of b.c. is "strong enough", 
the bulk free energies also yield information about normal fluctuations, as 
well as large deviations, in particle number and energy, inside large regions 
of macroscopic systems. The theory of such fluctuations is at the present 
time well developed for classical systems [0, 0, ^FS|, but almost nonexis- 



tent for quantum ones. It is the purpose of this note to make a beginning 
towards such a quantum theory. This is clearly desirable since the real world 
is quantum mechanical, with the classical description being an essentially 
uncontrolled approximation, albeit a very good one in many circumstances. 
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1.1 Classical systems 

We begin by considering a classical system of A^ particles of mass m in a 
domain, say a cubical box V C M'^, interacting with each other through a 
sufficiently rapidly decaying pair potential 0(r), e.g., a Lennard- Jones po- 
tential. The Hamiltonian of the system is then given by 

N N 

H{N,V;b) = —J2pl + ^ Yl 0M + $^^.(r.), (1.1) 

where pj G M*^, Tj G V , rij = |rj — fjl, and Mb(rj) represents the interaction of 
the i-th particle with the world outside of the boundary of V. This boundary 
interaction (indicated and in the sequel by b) is in addition to the action of 
the implicitly assumed "hard wall" which keeps the particles confined to 
V. The dynamic effect of the latter is to reflect the normal component of 
the particle's momentum when it hits the wall. However, sometimes it is 
convenient to replace it with periodic boundary conditions |p^'L| . 



For a macroscopic system in equilibrium at reciprocal temperature /? and 
chemical potential /i, the grand canonical Gibbs ensemble then gives the 
probability density for flnding exactly N particles inside V^ C M'^ at the 

phase point X^ = (ri, pi, . . . , r^r, p^v) = (Rat, Pn) G r^r = V^ x M*^^ as 

u{x^ I (3, /i, V, b) = mjwj) ' ^ ^^ 

where S is the grand canonical partition function 

oo 

= $^e^^^Q(/?,iV|y,6), (1.3) 

N=Q 

where Q(/3, N \V,b) is the canonical partition function. We use h'^^ , h being 
Planck's constant, as the unit of volume in the phase space Fjv, so A^ = 
h^ [3 / {2717X1) is the de Broglie wave length. The flnite- volume, boundary 
condition dependent, grand canonical pressure is 

p(/?,/i|\/,6) = /3-^|\/|-MogH(/3,/i|X^,6). (1.4) 
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Taking now the TL, V y R"^, we obtain, for a suitable class of b.c, 
an instrinsic (b.c. independent) grand canonical pressure p{(3,fi). This is 
related to the Helmholtz free energy density a(/3, p) obtained from the TL of 
the canonical ensemble, i.e., S is replaced by Q~^ in ( |1.4|) and the limit is 
taken in such a way that A^/|^| — > p, a specified particle density. The relation 
between p and a is given by the usual thermodynamic formula involving the 
Legendre transform 

p(/3, p) = sup[p/i - a(/?, p)] = 7r(/?, p), (1.5) 

p 

where 7r(/5, p) is the TL of the canonical pressure 

.(Ap) = -p« (1,6) 

and 

p(/3,p) = ^(/3,p) (L7) 

is the average density in the grand canonical ensemble. 

At a first order phase transition p(/?, p) is discontinuous and the left /right 
limits of the derivative on the r.h.s. of ( p..?] ) give the density in the coexisting 
phases. In our discussion we shall resctrict ourselves to values of the param- 
eters P and p where the system is in a unique phase. We can of course also 
go from the grand canonical pressure to the Helmholtz free energy density 
by the inverse of ( iLSf) , 

- a(/3, p) = sup[p(/3, p) -pp]. (1.8) 

Let P{Nv G A|l^| I /?, p, V, b) be the probability of finding a particle den- 
sity in V which lies in the interval A = [ni, 712], i.e., between the densities rii 
and n2. Then, for b in the right class of b.c, we have (almost by definition) 
that 

lim \V\'HogP{Nv e A\V\\ (3,p,V,b) = - mi[a{f3,n) - a{(3,p)], (1.9) 



where p is given by ( |1.7| ). In probabilistic language, this means that a(/5, n) is 
the "large deviation functional" for density fluctuations. (Note that a{(3,p) 
may be infinite for some values of p, i.e., when 0(r) = 00, for r < D, and p 
is above the close-packing density of balls with diameter D). 
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On the other hand, the fluctuations in all of V are clearly ensemble de- 
pendent (they are nonexistent in the canonical ensemble) and therefore not 
so physical. More relevant are the fluctuations not in the whole volume V 
but in a region A inside V. Of particular interest is the case when A is 
very large on the microscopic scale but still very small compared to V. The 
proper idealization of this situation is to first take the TL oi V /" M.'^ and 
then let A itself become very large. We are thus interested in the probability 
P{N\ G A|A| |/?,/i), for A a large region in an infinite system obtained by 
taking the TL of V. This probability should now be an intrinsic property of a 
uniform single-phase macroscopic system characterized either by a chemical 
potential /i or by a density p. 

A little thought shows that this probability corresponds to considering 
the grand canonical ensemble of a system of particles on a domain A with 
boundary interactions of the type 

oo 

M6(r,)=^0(|r,-Xfc|), r,Gy, XfcG^ (1.10) 

k=l 

i.e., we imagine that the boundary interactions come from particles of the 
same type as those inside A, specified to be at positions xi,X2, . . . outside 
A. These positions must then be averaged according to the infinite-volume 
Gibbs measure. It follows then, from the independence of the bulk properties 
of the system of the boundary conditions, that equation (|L9|) is still correct, 
that is 

lim |A|-MogP(iVAG A|A||/3,/i) = -inf[a(/3,n)-a(/5,p)]. (1.11) 



This relation is indeed a theorem for classical systems, under fairly general 
conditions g, ^, y. 

1.2 Quantum systems 

It is equation ( |1.11| ) and similar formulas for fluctuations in the energy den- 
sity which we want to generalize to quantum systems. To do this, we begin by 
considering the boundary conditions imposed on the A^-particle wave func- 
tions \E'(ri, . . . , rjv 1 1^) for a quantum system in the domain V. Usually this 
is done by requiring that whenever any Tj is at the boundary of V, Tj G dV, 
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then \1/ is equal to a times its normal derivative 

*(ri, . . . , r;v I ^) = « n,, ■ T^^(ri, . . . , r^ | y) (1.12) 

with a = corresponding to Dirichlet and a = +00 to Neumann boundary 
conditions. 

Denote by ba the elastic boundary condition ( |1.12| ). The existence of the 



TL of the grand canonical p{(3, /i | V, 60) has been proven for quantum systems 
with stable potentials 0, and for positive potentials it is established that 
the pressure does not depend on ba ||Ro|| . But, as far as we are aware, the 



dependence on Mfe(rj) has not been studied systematically, with the exception 
of the regime covered by the low-density expansion of Ginibre pi| , [BK|| . This 
only shows that the dependence on the boundary is not so well understood 
for continuous quantum systems. 

To investigate the density fluctuations in quantum systems we note that 
the momentum variables did not play any role in the derivation of ( |1.9|) 
and ( |1.11| ) for classical systems. The only thing relevant, when considering 
particle number fluctuations, is the probability density in the configuration 
space. This is given for a classical system by integrating u in ( |1.2| ) over 
the momentum variables, whose distribution is always a product of Gaus- 
sians (Maxwellians). For a quantum system, where the analog of ( |1.2| ) is the 
density matrix z/, the configuration probability density is given by the diag- 
onal elements of u in the position representation. For the grand canonical 
ensemble this can be written as 

W{K^\P,fi,V,ba) = ^"' '^ g ^' , (1.13) 

where \1/^ and E^ are the eigenstates and eigenvalues of Hn with the suitable 
statistics and ba b.c. [0, 0]. 

It is clear from the derivation of the TL |^, [T] that, when 0(r) is super- 
stable, the TL for the canonical ensemble exists for all p G [ni,n2] with b.c. 
ba- Then (|1.9| ) carries over to quantum systems. The real problem is how 
to prove ( p.. lip for these systems. W is no longer a Gibbs measure with a 
pair potential as interaction and there is no good reason to expect it to be a 
Gibbs measure for any other "reasonable" many-body potential |[vFS|| . (Even 



if the latter were the case, this potential would almost certainly depend on 
the density and temperature of the system and would therefore not carry 
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directly any information on ( |1.11|) .) It might in fact appear that there is 
no strong reason why ( p. . 1 1|) should hold for quantum systems. The reason 
for expecting it to be true is that it is a thermodynamic type relation and 
such relations are in general unaffected by the transition from the classical to 
the quantum formalism. More explicitly, we see the difference between ( |1.9|) 
and ( 1 . 1 1|) as involving only boundary type quantities which should become 



irrelevant when A is of macroscopic size. The proof of such a statement is 
however far from obvious (to us) and we therefore devote the rest of this note 
to proving it in the (technically) simplest case where there are no interactions 
between the particles, i.e., the ideal gas with either Bose-Einstein or Fermi- 
Dirac statistics. It turns out that even in this case the proof requires a certain 
amount of work. 

To finish this introduction, we note that the same reasoning which leads 
to ( |1.11| ) also gives the well known result that the variance of N\, divided 
by |A|, is given, for |A| y oo, by the compressibility ((A^a — p|A|)^) /|A| —> 
(3~^{d'^p/diJ,'^){l3, fi). Furthermore, a central limit theorem for the random 
variable ^ = lim|A|/'oo(^A — p|A|)/a/p|A| holds. These results are also ex- 
pected to remain valid for quantum systems and are proven here in the 
non-interacting case. 



2 Main results 

We consider a d-dimensional square box V = [—i/2, Ijlf . For computational 
convenience we choose periodic boundary conditions. The infinite-volume 
thermal state we construct below does not depend on this particular choice 
of the boundary conditions, however ||BK|| . In V there is an ideal fluid (either 
Fermi or Bose) in thermal equilibrium, as described by the grand canonical 
ensemble. We label the Bose fluid, shorthand BE, with the index + and the 
Fermi fluid, shorthand FD, with the index — and introduce the Fock space 

oo 

^I = C©0L^(n, (2.1) 

n=l 

where L^(^") is the ra-particle space of all symmetric, resp. antisymmetric, 
square-integrable functions on V'^ . Of course, for ra = 1, L\{y') = L'^{V). 
In the sequel, in order to keep the notation light, we will often drop sub- or 
superscripts whenever there is no ambiguity. 
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Particles do not interact. Therefore the many-particle Hamiltonian in the 
box V can be written conveniently in the form 

oo n 

^^ = G^y^ J^-'-^hy®---^! , (2.2) 

~ ~ «*'' position out of n 

where hy, the one-particle Hamiltonian on L'^{V), is defined through the 
one-particle energy e{k) in momentum space. This means that, if \k) denotes 
the momentum eigenvector (represented in L'^iV) as ^y (x) = e*'^'^), then 
hv\k) = e{k)\k) with keV = (27rZ/£)"', the dual of V. 

We assume e{k) to be continuous, e(0) = as a normalization, and e{k) > 
for fc 7^ 0. Also e{k) ^ \k\'^ for small k and e{k) > \k\°' for large k, with 
a, 7 > 0. Furthermore we require 



I^'^I 



d^ke'''-'' ^ 



/3e(fc)-/3M _ 



< OO (2.3) 



ioY e = ±1, P > 0, and suitable fi. 

The standard example of a non-relativistic, resp. relativistic, kinetic en- 
ergy for a particle of mass m is e{k) = k'^/{2m), resp. e{k) = \Jm^c'^ + k'^(? — 
mc? (having set Planck's constant h= \\ Both functions satisfy the above 
conditions. The relativistic case includes m = 0, although this is not imme- 



diately obvious, cf. Appendix |A.1| for details. 

We observe that Hy may be rewritten as a quadratic form in the creation 
and annihilation operators on the Fock space 5". Let a\. be the operator that 
creates a particle in the state \k) and a^ the corresponding annihilator. Then 

Ey = ^e{k)alak = ^{j\hy\k)a*ak = {a\hy\a) . (2.4) 

k j,k 

We fix /? > and yU G M for FD, resp., ;U < for BE. The grand canonical 
state in the volume V is defined by 



,v 



Tt^v {a e-f^^^+f^^^) 



(^)L = " ^v... (2-5) 






if^) 



for every bounded operator A on 3^^. N = Ny is the operator for the 

number of particles in the box V, N\i2 f^yn^ = nl/,2 (y„), and S^(/i) = 
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Trgrv(e~^^^^^^^) denotes the partition function. As is well known (see, for 
example, |0) we have 

Elifi) = 1[{1 - e-f^^^+f'^y^ (2.6) 

k 

Hr(/i) = IKl + e-f^''^-^'"^) . (2.7) 

k 

The infinite- volume thermal state is defined through the limit 

(■)= lim (■)^, (2.8) 



when taking averages of local observables PR| , Sec. 2.6]. 

Taking the infinite volume limit of ( |2.6| ) and ( |2.7] ) one obtains the grand 
canonical pressure 

p^(^) = lim i^lP^ = ---^ / d^k log (1 - 6 e-^^(^)+^^) (2.9) 

and the average density 

p(^) = '^u\ = _}_ [ d^k , ,,, \ . (2.10) 

p_ is real analytic on the whole axis, whereas p+ is real analytic only for fi < 

and has a finite limit as /i — > 0_. For convenience, we define p+{n) = oo for 

/i > 0. The slope of p+ at 0„ is related to the Bose-Einstein condensation. 

We set 

1 /" 1 

p^ = p(0_) = - — (fk-—. . (2.11) 

'^ ^^ ' {2-kY J eP<^)-l ^ ' 

By the properties of e(/c), pc = oo ior d < 7, and is finite otherwise, pc is the 
maximal density of the normal fluid and any surplus density is condensed 
into the /c = ground state. To simplify the notation we use pc also in the 
case of an ideal Fermi fluid, setting it equal to 00. 

The infinite system is assumed to be in a pure thermal state, obtained 
through the limit (|2.8|) at the reference chemical potential p. In this state 
the average density is p = p{p) < pc- We define the translated pressure by 

^,,^(A) = QeiX) = Peip + A) - peip) . (2.12) 
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Qs is convex up, increasing, g{0) = 0, and g'^{0) = p. For large negative values 
we have 

lim 9ei\) = -Peif^) , lim ^^(A) = 0, (2.13) 

A— >— oo A^— oo 

whereas for positive values 

lim g_{\) = cx) , lim g'_{\) = oo (2.14) 

A^oo A^oo 

in the case of fermions and 

lim ^+(A) = p+(0) - p+(/i) , lim ^^(A) = pc (2.15) 

A— >— /I A— >— /I 

for bosons, with g+{\) = oo for A > —p. 

We define the rate function /^ as the Legendre transform of gs, i.e., 

fsA^) = fe{x) = inf (fl'e(A) - \x) = ge{\o) - Kx. (2.16) 



Here Aq = Xo{x) is the minimizer of g{X) — \x, which is unique by convexity. 
For a; < we have Aq = — oo. For < x < p^, it is determined by g'{\o) = x, 
while for x > pc we have Aq = —p. This shows that f{x) = — oo on the 
half-line {x < 0} and f{x) finite elsewhere. In particular, / is convex down, 
strictly convex for < x < pc, and f+{x) = p{0) — p{p) + px, for x > pc, as 
a trace of the Bose-Einstein condensation. 

Let us now consider a small subvolume A of our (already infinite) con- 
tainer V. The precise shape of A plays no role, only the "surface area" 
should be small compared to its volume |A|. Thus, by A /" M^ we mean a 
sequence of subdomains such that for each A there exists a subset A' of A 
with |A'|/|A| -^ 1 and dist(A', M^ \ A) ^ oo. 

Let A^A be the number operator for the particles in A. With respect to 
( ■ ), A^A has some probability distribution. We follow the usual practice and 
use the same symbol A^a to denote also the corresponding random variable. 
Its distribution is indicated by P, averages again by ( ■ ). 

We are now in a position to state the main result. 

Theorem 2.1 Let P > and p < for BE, resp. p E M. for FD. Consider 
an interval I = [a, b] . Then we have, for a < Pc, 

hm -^logP^A^A G |A|/}) = sup4^(a:) , 
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and, for a > pc, 

limsup-^logP({iVA e |A|/}) < sup/,,^(a;). 



3 Large deviations in the density 



In this section we explain how Theorem p.l| follows from the asymptotic 



behavior of the generating function (^e^-^^^ ^ . 

Lemma 3.1 There exists a \max{^) such that i^e^^^'^^ < oo for all A < 

Amaa;(A) and {^e^^^^)^ = oo for all A > \max{,J^)- For FD we have Xmax{J^) = 
cxD, whereas for BE Xmaxi-^) < oo with Amax(A) \ — /i, as A y M.^. 



Theorem 3.2 The limit 



lie. )3|A| ^' =9.W, (3-17) 



including any finite number of derivatives, exists uniformly on compacts of 
M for FD, resp. of (— oo, — /x) for BE. 

These results are proved in Sections ^ and |[ 

Inferring Theorem ^TI| from our information on the generating function is 
a standard argument from the theory of large deviations 0, |0l . 
The probability of the event in question can be rewritten as 

QA=(X|A|/(iVA)>, (3.18) 

where xa is the indicator function of the set ACM. To make this event 
typical we introduce the modified average 

{■), = ^{-e'''''), (3.19) 

^A 

where A < Xmax and the partition function Zx = (e^'^^^). With respect to 
this new state, (|3.18| ) can be expressed as 

QA = Zx{e-f''''^X\A\iiNA))^. (3.20) 
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An upper bound for Q\ comes from the exponential Chebychev inequality, 

Qa < (^e^MNA-a\A\)^ ^ z, e-^^'^l^l (3.21) 

for any A < Xmax- 

The lower bound requires more effort and can be carried out provided 
a < pc- One uses ( p.l3|) -( pTT5D to show that there exists a Aq < Xmax such 
that g'{Xo) = ct- Differentiating (|3.17|) twice w.r.t. A we obtain 



lim 



(Na) 



Ao 



9'{K) = p{p + Xo) = a, 



lim — - 

jd |A| 



|A| 

(iVl)A.-((A^A)J 



1^^^+^°^ 



(3.22) 
(3.23) 



which is finite. This means that the event {N/<^ ^ |A|a} is typical for the new 
state and a law of large numbers holds. Notice that Aq > 0, since p is strictly 
increasing in p. From (|3.20|) , Vc G (a, b), 



Qa = ZA„(e-^^°^-X|A|[a,c](^A)>,^ 

> Z,„e-'^^-l^l(x|A|M(A^A)>,„ 



(3.24) 



for some a E (0, 1) and |A| large. In fact, (x|A|[a,c](A^A));^ -^ 1/2 as A /^ M°' 
Therefore, when a < pc, v^e obtain from ( |3.21|) , ( |3.24| ) and Theorem ^ 



^(A,) - Xoc + o(l) < i^^ < giXo) - Xoa + o(l) 
Since c G (a, h) is arbitrary, we conclude that 

lim —^TTT- = g{Xo) - Xott = f{a) = sup f{x) , 



(3.25) 



A/'K'* f3\A\ 



xS:[a,b] 



(3.26) 



where / is the rate function defined in (|2.16|) . Ao is the same as in the defi- 
nition of the Legendre transform ( p.l6[ ), because of ( p.22| ). The last equality 
comes from the convexity of /. 

The second assertion of the theorem is now easy, since the right inequality 
in (|3.25| ) holds without any restriction in a. Q.E.D. 

Theorem p.l| implies the central limit theorem for the density in A. 
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Corollary 3.3 Under the assumptions of Theorem 2A, the moments of the 
variable ^A = (^a — (^A))/|Ap^ converge, as A y W^, to those of a Gaussian 
with variance (3~^{dp/d^){^). 



Proof. The fc-th cumulant of ^a is given by 



CA(fc) 



1 



^fc|^|fc/2 



d 



k 



dX^ 



\og{eP^^^) 



(3.27) 



A=0 



k>2. From Theorem U, Ca{2) -^ /9-^^"(0) = (3-\dp/dn){fi), whereas, 
for k > 2, C^ik) — >■ 0. Also Ca(1) = 0. These limits are the cumulants of a 
centered Gaussian variable with the specified variance. Q.E.D. 



4 Generating function 

We derive a determinant formula for the generating function {^e'^^^"^) . With 
its help we prove the claims of Lemma \i.l[ We will see in the next section 
that is convenient to introduce the variables ( = e^'^ and ( = ( — 1. 
By dU), we have 

- (3 Hv + (3 pNv = {ai\{- (3hv + f3ijlv)\ai) = (a|Av|a), (4.1) 

/3AiVA = (a|/3AxA|a) = (a|5A|a), (4.2) 

which implicitly define Ay and B^ as linear operators on L'^{V). We will use 
the following identity. 

Lemma 4.1 Let A,B be self-adjoint and bounded from above. Then there 
exists a self-adjoint operator C such that e^e^e^ = e*" and 

g(a|A|a>g(a|iJ|a>g(a|A|a> ^ g(a|C|a) 

for both BE and FD. 

Proof. See Appendix. 

We apply Lemma ^?T| with A = Av/2 and B = B^, after a symmetrization 
of the density matrix in ( |2.5|) . Then, using also definition ( p.8|) , 

(e"-") = hm :5jd£^l^ (4 3) 

^ ' y/RdTr^^(e<'^l^^'l^)) ' ^ ' ^ 
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Evaluating the trace of a quadratic form in a*, aj is a standard calculation for 
both BE and FD. Let us consider first the case of fermions. For a self-adjoint 
operator A on L'^iV) such that e^ is trace-class, we have 



Tr^v(e<"l^l">)=dety(ly + 



det(l + xye^Xv), 



(4.4) 



where dety is the determinant on L'^{V) and det the determinant on L'^{W^). 
Here and in the sequel we refer to the theory of infinite determinants, as 
found, e.g., in ||RS| , Sec. XIII. 17]. e^^ is obviously trace-class, and so is e*", 
since e^^ is bounded. Using the definition of C, we obtain 



dety(lv + e^) 
dety(ly + e^^; 



Av/2Ra^Av/2 



dety [(ly + e^^)-^(ly + e^^/^e^'^e 

dety [ly + (ly + e^^y'e^^^^e''^ - ly)e^^/2] 



det 



1 + C Xa ^y- Xa 



(4.5) 



where Dy. = (l + e^^)-^e^^. We used the fact that e^'^ = (e^^-l)xA + l = 
Cxa + 1 and the cyclicity of the trace in the definition of the determinant. 
Finally, from (^) and ( ^31) , 



(e'^^^A), 



lim det 



1 + C XA ^y,- XA 



(4.6) 



One would like to take the limit on V inside the determinant by replacing 
Dy- with the corresponding operator on L^(R'^) defined as 

{D^){k) = Z{k)ip{k), {D^ij){x) = jdyd^y- x) ^{y), (4.7) 

where ^ denotes the Fourier transform and 

1 



d_{k) 



1 + (^P{e{k)-ti) ■ 



(4.8) 



Notice that (i_ G L^{M.'^) by our assumptions on e{k) and so d- G L°°(M'^). 
Moreover, ( |2.3| ) ensures that d^ G L^(M'^). 

By |[RIJ| , Sec. XIII. 17, Lemma 4(d)] one has to estabhsh that XA-Dy,- Xa 
tends to xa D_ xa in the trace norm. 
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Lemma 4.2 Let d be a continuous integrahle function on M.^ . We define 
D through ([4.7|) as a linear operator acting on L'^iW^). Furthermore we 
define Dy by Dv\k) = d{k)\k) on L'^{V) and by Dy = on the orthogonal 
complement L'^iW'' \ V). Then, for A C V^, Xa-^vXa '^''^d Xa-Dxa o^re trace- 
class, and 

lim Tr|xA(-Dy - D)xa\ = 0. 



Proof. See Appendix. 
We conclude that 

(r->_ = det(l + CXA/^-XA), (4.9) 

with ( = (- 1. 

For bosons we proceed in the same way, except that ( [4.4|) is replaced by 

Tr^v(e<-l^l->) = dety(ly-e-^)-\ (4.10) 

requiring in addition ||e"^|| < 1. In fact, for ||e^|| > 1, the l.h.s. of (|4.1(J|) is oo. 



whereas the r.h.s. might be finite if 1 is not an eigenvalue of the trace-class 
operator e^. In our case, by assumption ||e^^|| < 1. As for e*-^, the function 

A ^ lle^ll = II (e^^ - l)e^^/2xAe^^/' + e^^|| (4.11) 

is increasing and Xmaxi-^) is defined to be that A which makes it equal to 1. 
Since the r.h.s. of ([4.11|) is increasing in A and its sup is e^'^||e'^^|| = e^*-'^"'"^\ 



then one checks that Xmaxi-^) \ - 
computation for FD, we have 


-fi, as A / W^. 


Therefore, following the 


(C^) - lim det(l-e'^) 

^ /+ y/Rddet(l-e^v-) 


-1 
-1 


det(l + CxA/^+XA)"', 


(4.12) 


for A < Xmax and cxd otherwise. Here D+ 
is defined as in (4.7]) with 


, the limit of Dv,+ = (e^^ - 


-ly'e^y, 


d;{k) = 




1 




(4.13) 


1- 


g/3(.(fc)-M) • 


Equation (4.12|) is the analogue of 


([4.9) 


and proves 


Lemma 3.1. 
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5 Infinite volume limit 

Instead of the chemical potential, in this section we use the fugacity z = e^'^, 
regarding it as a complex variable. This will come out handy for the proof 
of Theorem |3.2| . The variables ( and (, defined at the beginning of the 
previous section, will also be extended to the complex plane. In this setup 
the translated pressure ( |2.12| ) becomes 



9ziC)=PizC)-piz), (5.1) 

where, with a slight abuse of notation, we keep the same name for the pressure 
as a function of the fugacity. 

Expressions (|2.9| )-( prT0| ) for the pressure and the average density define 
two analytic functions of fi in 

E+ = {Re fi < 0} U {Re fi > 0, Imfiy^ 271 J / (3, \f J eZ}, (5.2) 

E_ = {Re/i<0}U{Re/i>0, Im/i^ (2j + l)7r//?,Vj gZ}. (5.3) 

Hence QeiC) is analytic in 

G+ = C\[z-\+oo); G- = C\{-oo,-z-^]. (5.4) 

We proceed to give the proof of Theorem ^.2[ Let K G G^ he a. compact 



set in the complex plane. We choose K such that L = K f] M"^ is also 
compact, since its image through the function C '"^ -^ verifies the hypotheses 
of the theorem. Our argument, however, is valid for any K. Without loss of 
generality, we can assume that 1 E K. 

For ( restricted to Gs fl M"*", let us define 

0t(C) = m 1°S (^"""^s . = -m^^ 1°S(1 + CxaD^xa) (5.5) 

A ^'^ A 



according to ( [4.9| ) and ( [4.12|) . The proof of Theorem |3.2| will be subdivided 



into three steps. 

1. (j)^ can be analytically continued to G^. 

2. There is a positive r such that 0^(C) converges uniformly to PgziC) for 
|C-1| <r. 



Large deviations for ideal quantum systems 17 

3. \(j)^\ is uniformly bounded on K. Therefore by Vitali's lemma |T], 
Sec. 5.21] 10^1 and any finite number of its derivatives converge uni- 
formly on K. 

Step 1. We leave the proof of the following lemma for the Appendix. 

Lemma 5.1 The function 0^(C)j os defined by the tri 
in Gs- 

Step 2. Expanding the log in (|5.5| ) one has, for \(\ < \\D^^^^^ 



ell ) 



oo 



m— 1 



0f (C + 1) = -^Tr J2 ^-^{CxkD^XkT. (5.6) 

|A| ^ m 



We would like to interchange the summation with the trace. To do so, we 
need dominated convergence for the series: 

icxAi^.xAr < icrp.ir-'(xA^.XA). (5.7) 

Since |A|^^Tr(xA-D£XA) = d^{0) (see proof of Lemma ^^ in the Appendix), 
each term of ( p. 61 ) is bounded by a term of an integrable series independent 
of A. Therefore, for the same C's as above, 

</>^(C + 1) = -^ E -^ — ^ m ^^XADsXAr. (5.8) 

m=l ' ' 

Suppose that we are able to prove that 

hm -^ Ti{xAD,XAr = [ dk[lik)r, (5.9) 

A/Rt* |A| J 



with a rest bounded above by mR"^ for some positive constant R. Then, 
using ( |4.8| ) and ( [4.13[ ), we would have that, for any r < min{||De||^^, i?~^}. 



uniformly for \(\ < r, 

^ °° f 1 \in—l/-m r / 1 \ ' 

lim0^(C + l) = -sY^-^ ^— dki- ---^ 



771=1 
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the last equality coming from (|Z!9|). This would complete Step 2. 
Let us pursue this project. One sees that 

/ dk[de{k)]"' = {d, * d, *■ ■ ■ * d,l {0) (5.11) 

m times 

= -— dxi dx2ds{xi-X2)--- dx^d^{xni-^i-Xm)d,{x^-xi). 

The normalized integration over xi is harmless since, by translation invari- 
ance, the integrand does not depend on that variable. On the other hand, it 
is not hard to verify that 

±-TrixADeXAr = m/^^i {xi\iXAD,XAr\^^i) (5.12) 

= 777 / '^^1 / dX2ds{Xi-X2)--- / dXmde{Xni-l - Xm) de{Xm - Xi). 
W J A J A J A 



In view of ( |5.9| ), we want to compare (|5.11| ) with (|5.12|) . We observe that 




A J A J A 



m—l 



i=l 




E I- 1 IL- b (513) 



m— 1 times m—l times i— 1 times m—l—i times 

Subtracting (|5.11| ) from ( p.l2| ) leads then to m — 1 terms of the form 
1 

'A JA' 



, dXi I dX2 de{Xi - X2) ■ ■ ■ I dXmde{Xm-l-X.m)de{Xm-Xi), (5.14) 



where the sets A^, . . . , Am can be either M°' or A. ( ^.14| ) holds because, due to 
the cyclicity of the integration variables, one can cyclically permute the order 
of integration without touching the integrand. We overestimate by switching 
to absolute values and integrating x^, . . . ,Xm over M.'^, 

-— / dxi {xA--xAd,\) * 141 * ■ ■ ■ * |4|)(xi) = T-TT dxiUA{xi), (5.15) 

W J A VM J A 

which defines u^^xi). To estimate this function, we use recursively the rela- 
tion ||/*5'||oo < ll/lloollfi'lli and obtain 

ua{xi) < IMeII™""^ sup 141- (5.16) 

A^-xi 
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Recalling now the definition of A' given before the statement of Theorem pTT 



one sees that, if xi G A' and y G A^ — xi, then \y\ —>■ oo as A y M . Hence, 

from ( |5.16| ), 

sup sup ||(ie||]^'"^"'^'UA(xi) ^ 0. (5-17) 

xiGA' m>l 

Also from ( 5.16| ), pointwise in xi, 

sup ||4||];'"+^nA(a;i) < ||4||oo- 

m>l 



(5.1^ 



When we average over xi G A, the last two relations and the properties of A' 
prove that 



lim sup nielli 



m>l 



lAI 



dxiUAixi) = 0. 



(5.19) 



This takes care of each term as in ( p.l4| ), and we have m — 1 of these terms. 
Hence (|5.9| ) holds with R= \\d^\\i- This ends Step 2. 



Step 3. Again we expand (|5.5| ) in powers of (, but this time about a generic 
Co^Ge-l (see (U)). We obtain 



i-Tr log(l + CXaDsXa) = ^Tr log(l + CoXa^.Xa) (5.20) 

+ m^^E^^ Ul + CoXADsXAr'XADeXA) (C " Co)"^- 

' ' m=l 



Let us estimate this series. First of all, using some spectral theory ||W| , 
Sec. 7.4], 



< 
< 



dist(l,cr(-CoXA^eXA)) 
-1 



:1 + CoXaD,xa)-' 

distil, cr{-CoDe)) 
since we know from definitions (|4.7| ), (|4.8| ) and (|4.13|) that 



(5.21) 



aixAD^XA) C [0,||xa/^-Xa||]C [0,||D_||] 

= aiD_) = [0,l/{l + z-% 

aixAD+XA) C [-||xa/^+Xa||,0]c [-||D+||,0] 

= aiD^) = [l/il-z-'),0]. 



(5.22) 
(5.23) 
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Repeating the same reasoning as in Step 2, we use the above to exchange the 
trace with the summation in ( |5.20|) — which is legal for small |C ~ Col as to be 
determined shortly. This yields a new series, whose m-th term is bounded 
above by 

dist(i, a(-Co/^,))j iiD,ir-M,(o) ic - Cor = a (6(Co) ic - Coi) , 

(5.24) 
where ^^(O) = |A| ^"Yiixt^DeXiv)- Hence, in view of (p.5|), ( [5.20| ) implies 

I0^(C + 1)1 < l0^(Co + 1)1 + a-^%fc% < |0^(Co + 1)1 + a, (5.25) 

l-KCo)|C-Co| 

for|C-Co|<(2KCo))^'. 

The crucial fact is that &(C)"^ stays away from zero when C is away 
from the boundary of G^ — 1. This can be seen via the following argument, 
exploiting ( [O^ and (|^-(|]23l). In the FD case a(-Co£'-) is a segment 



that has one endpoint at the origin and the phase of —Co is the angle it 
forms with the positive semi- axis. This means that, as long as Co does not 
go anywhere near the negative semi-axis, we are safe. For Co ^ (~^~^~ 1? 0) 
(see ( ^.4]) ), cr(— Co-D_) is contained in IR+. However, notice from ( |5.22| ) that 



the other endpoint is located at — Co/(l + 2^^) < 1- For BE the reasoning 
is analogous, except that in this case the phase of Co is the angle between 
(j(— C0-D+) and M^. Therefore the "safe" span is the complement of the 
positive semi-axis. Also, if Co ^ (0,-2"^— 1) (again see (|5.4])), the "floating" 
endpoint of cr(— ^o-D+) is found at (^qI{z^^— 1) < 1. 

With the above estimate we can use ( ^.25| ) recursively. If |Co| < ?", from 
Step 2, |0^(Co + 1)1 < ^5 for some M, since 0^ converges uniformly there. 
Then, from ( |5.25| ), we have that |0^(Ci + 1)1 < M -|- a, for any Ci such that 
ICi - Col <_(2&(Co))"^- Proceeding, we see that |0^(Cfc + 1)| < M + ka, 
whenever |Cfc — Cfc-i| < (2&(Cfc-i))~^- Iii this way we will cover K in finitely 
many steps^since it keeps at a certain distance from the boundary of G^ 
and the {b{C,k))~^ are bounded below. This completes Step 3, i.e., 0^(C) is 
bounded on K and Vitali's lemma can be applied. Q.E.D. 
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A Appendices 

A.l Relativistic massless particles 

We prove that the energy dispersion e{k) = c\k\ satisfies our assumptions. 
The only condition to be checked is (|2.3|) , that is, the Fourier transform of 
k H- »• (^e^^^\'^\~^^^ — e)~^ is in L^(M'^). This is a consequence of the following 

Lemma A.l Let f : [0, +00) — > C be of Schwartz class. With the common 
abuse of notation, denote by /(|^|) the Fourier transform of f{\x\), for x,^ G 
W^ . Then, for some positive C , 



fm)< 



c 



\i\ 



d+l' 



Proof. For simplicity let us write ^ = |.^|. The Fourier transform of a radial 
function is 

m) = ^-Tjlir / drfir)r'/'j,;,_M) , (A.l) 



^'/'-' JO 
cf. ||SW| , Chap. IV, Th. 3.3], where J^ is the standard Bessel function of order 



u [|Wa|| . One has 



X 



^"<^' " ¥tU7TY,- (^-2' 



for a; ^ 0, whereas 

Mx) = \ — 
V irx 



cos ( X n J + g[x) 



(A.3) 



with g{x) — > for x ^ 00. Using the relation 

dtfJ^.i{t)=x''J^{x) (A.4) 



we integrate ( |A.1| ) by parts repeatedly, taking into account also ( |A.2| ) and 
the hypothesis on /. After n integrations we get, up to constants, n terms 
of the form 

___y^ rfr/»(r)r'^/2-"+V,/2+„-i(rO, (A5) 

with i = 1, . . . , n. For our purposes it suffices to iterate up to ra > d/2 + 2. 
In fact, if i is such that d/2 — n + i > —d, then in (|A.5|) we can estimate the 



22 J. L. Lebowitz, M. Lend and H. Spohn 



Bessel function by a constant. The integral converges by the rapid decay of 
f^^' and the whole term is of the order ^^'^^^ or better. For smaller values of 
i, the estimate uses ( |A.2| ), for x G [0, a], and ( |A.3D otherwise. Since \f^^^\ < c, 
( |A.5| ) is bounded by 







fd/2+n-l 



the second integral being convergent because of the choice of i. Q.E.D. 



A. 2 Proof of Lemma 4.1 

As before, we set e = ±1, according to either bosons or fermions. A general 
/ G L^iV) can be expanded in the Fourier basis as / = ^^/fclfc)- The 
corresponding creation operator is then defined by 

fcev 

For the sake of simplicity, we denote A = (a|y4|a) = J2ij ^ij^t^j (same for 
!B). Recalling the canonical (anti) commutation relations. 



one calculates that 
and in exponential form 



ai,a*]_e = aia*-ea*ai 


= Sij] 


ttiyttj S = ttittj — e ttjtti 


= 0, 


[A,aifr] = aiAf)\ 





(A.8) 
(A.9) 

e*-^a(/)*e-*^ = a(e*^/)*. (A.IO) 



Now, let |0) be the ground state of 5"^. For n E N, and /i, f2, ■ ■ ■ , fn G 
L'^{V), the finite linear combinations of the states 

l/i, /2, . . . , fn) = aihYaihT ■ ■ ■ a(/„)lO) (A.ll) 

are dense in 3", which is another way of stating that |0) is cyclic w.r.t. the 
algebra generated by the creation operators. Therefore, we need only test 
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our assertion on vectors of the type ( A.ll| ). Using ( A.IC ) with t = 1, and 
observing that yi|0) = 0, we obtain 

e^l/i,...,/.) = eMfiTe^-^-'-eMfnTe-^lO) 

= a(e^/i)*---a(e^/„)*|0) (A.12) 

= |e^/i,...,e^/„). 

The existence of C is a consequence of the spectral theorem. We call C the 
corresponding quadratic form in a*, aj. Through the repeated use of ( [A.12| ), 
one checks that applying c'^e'^c'^ to the states ( [A.ll|) is the same as applying 
e'^. The semiboundedness of A and B ensures that the domain of their 
exponentials is the whole L'^{V) and all quantities are well defined. Q.E.D. 

A. 3 Proof of Lemma 4.2 

For any symmetric operator A, xaAxa < XaI^Ixa- Hence |xa^Xa| < 
XaI^Ixa and Tr|xA^XA| < TrA|A|. When A = Dy, the convergence of 
the trace is proven by writing the further estimate TrA|Dy| < TrylD^I 
and then summing an intregrable sequence of discrete eigenvalues. For 
A = D, one uses the Dirac-delta representation of the trace to find out 
that TrA|-D| = |A| (27r)^°' J |c?|. The first assertion of the lemma has been 
proven. 

As for the second part, let us write 

Tt{\xa{Dv - D)xa\) = TT{UxA{Dy - D)xa) (A.13) 

= Try(f/xA^v'XA) - Tr{UxADxA) =Ti>-T, 

where U is the partial isometry -^^(A) — > L'^i^) that realizes the spectral 
decomposition as in |[RS| , Th. IV. 10]. It is convenient to use the position 
representation for the bases. So, ip^''\x) = e^^'^ and, as defined in Section]^, 
V'y = V'^^''xy- Let us work on Tf. using the cyclicity of the trace one obtains 



T. = Y,Y.^i^'i\xKUxADy\i^^P) 
kev 

= llT.^ik){i^^v^\XAUxA\i^^v) (A.14) 

kev 

= llT.^ik){^P^'^\XAUxA\^P^'^), 



kev 
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the last equality being due to the presence of the indicator functions xa- In 
complete analogy with the above, 

T = j^^ J dkd{k) Oj^'^IxaUxaH^^'^). (A.15) 



Since \{iP^''^\xaUxa\¥'''')\ < |A|, it is obvious that i^^J§ tends to J ^KT^ ) for 



oo. Q.E.D. 



A. 4 Proof of Lemma 5.1 

With regard to ([4.9|) and ( |4.12|) , det(l + (xADeXA) is entire in ( (hence in () 
by ||RS| , Sec. XIII. 17, Lemma 4(c)]. In order to evaluate its log (on the suitable 



Riemann surface) we need to avoid the zeros. Using ||RS| , Th. XIII. 106], we 
want to make sure that cr {—CxaD^xa) does not hit 1. Step 3 in Section |^ 
(see in particular formulas ( |5.22| )-( p.23[ ) and the last paragraphs) shows that 
this is never the case if C ^ (~oo, —z^^ — 1], for FD, or ^ ^ [z~^ — 1, +cxo), 
for BE. Q.E.D. 

Actually, we can say more. Consider FD, just to fix the ideas. We see from 
(|5.22| ) that the "floating" endpoint of cr(— Cxa-D-Xa) is strictly contained in 
the segment (0, —(^/{l+z~^)), which means that ( is allowed to exceed slightly 
G- — 1, as given by ( |5.4| ), without any vanishing of (^]9|). The above, and 
an analogous argument for BE, prove that for each flnite A the domain of 
analyticity of 4>^{C) is indeed strictly bigger than G^- 



For the bosonic case this fact is related to Lemma |3.1| . In a few words, 
XaD+Xa can be thought of as deflning a Hamiltonian /i^ in L^(A), via the 
relation 

XaD+xa = (e-^'^A+^^^^ - 1)-^ e-'^'^A+/3MiA_ (a.16) 



(Compare this with the deflnition of -0^,+ given after formula ( |4.12| )). Then 



all the calculations we have carried out in Section ^ are about the grand 
canonical ensemble of a system of bosons on A, with energy operator /i^ 
and chemical potential A + /x. The ground state of /i^ is strictly positive, 
in analogy to the Hamiltonian with Dirichlet b.c. In that case the upper 
bound for the chemical potential, Amax(A) + /i, is stricly bigger than zero 
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A. 5 The case A = V 

We discuss the less physical case of the large deviations for the density in the 
"large box" V. 

Instead of first taking the limit \^ ^ R'^ of the state ( ■ )„ , and then 
looking at the asymptotic properties of the quantity A^a/|A|, for A y W^, 
we now consider the large deviations for the variable A'^v'/I^l w.r.t. (■) , 
subject to the single limit V y M.'^. We call this the case A = V. 



We immediately see that Theorem ^^ is a trivial identity in this setup 



l3XNv\^ 



lim 



log(e^^^^) 

m 



lim 



logS^(/x + A)-logS^(/i) 

/3\V\ 



(A.17) 



Therefore the equivalent of Theorem |2.1| follows in the same way as outlined 
in Section |. 

In the case A = V, however, it turns out that we can do more than just 
this. At least for some interesting cases, we can provide a lower bound for the 
large deviation relation even in the BE condensation regime, i.e. for a > pc- 
This is how it is done. 



With an eye to Section H, we see that, due to (|2.15| ), the problem is that 
there is no fixed Aq that verifies ( |3.22| ). In other words, no fixed value of the 
"extra" chemical potential A can be found such that the modified state 



(■); 






(■) 



fi+X 



(A.18) 



has an average density a (compare the above definition with ( p.l9| )). On the 
other hand, as it is customary in the theory of Bose-Einstein condensation 
(e.g., | vLP| and references therein), one can take a variable Xy such that 



{Nvh 
\V\ 



p^{ix + \x 



d logH^(/i + Ay) 
dii (3\V\ 



a. 



(A.19) 



This formula is the analogue of ( p. 221) . As regards the analogue of (|3.23|) , it 
is not too hard to see that 



\v\ 



(^V>A. - i^Ny)J 



dp 
dp 



V 



■(/^ + Ai 



IV^I. 



(A.20) 
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That is, the variance of A'"y/|\^| is of the order of a constant, at the hmit. 
This means that our new state makes the density a average but not typical, 
in the sense that no law of large numbers holds. Thus, we cannot apply 
estimate (|3.24|) tout court. 

However, we do not really need the density a to be typical. The afore- 
mentioned argument goes through all the same provided that, for any c < a, 
we find an a G (0, 1) such that {x\v\[a,c\{Ny)'^ ^ > a, for V big enough. 

This suggests that we look at the asymptotic distribution of iVy/ |V^| in 
the grand canonical ensemble, when the infinite-volume on V is taken under 
the restriction (|A.19|) . This is a known object for several choices of V, e{k) 



and d. For instance, when \^ is a three-dimensional box (periodic, Neumann 
and Dirichlet b.c. apply) which expands isotropically and the energy is the 
non-relativistic e{k) = k'^/{2m), it is called the Kac distribution. Details of 



the computation can be found in [ZUK|, whereas accounts or generalizations 



of the result are presented in many other works [j^, ^EP, vLL]. It turns out 
that, for a > pc, 

lim (xwiANv)). = -^— [ exp (-^^^1 dx. (A.21) 

Therefore the above claim holds true under certain assumptions and, in anal- 
ogy with ( p.25| ), we can write that, for V y W^, 



with Ay —>■ — yU_ and c arbitrarily close to a. Hence 

^.^ logP({iVye^^^ > g(-p)+^,a = /(a) = sup fix). (A.23) 

y/Rd fi\V\ ^(.[afi] 

We have thus shown that in the case A = V, the rate function for the 
large deviations in the density exists and is equal to the Legendre transform 
of the translated pressure, even in the BE condensation regime. 
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